Phase space analysis for three and four massive particles in final states 
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We propose formulae for computing the phase space integrals of 1 — > 3 and 1^-4 processes with 
massive particles in final states. As an application of these formulae we study the final state mass 
effects in some interesting phenomenological cases, giving fully integrated analytic results for the 
corresponding phase spaces. We consider also the Bs — Bs process at NNLO and calculate one of 
the most complicated master integrals, which contributes to the AFb^ at 0{al). 
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I. INTRODUCTION 

The Bs,d -meson rare decays and oscillations are known 
to be a unique source of indirect information about 
physics at scales of several hundreds GeV. In the Stan- 
dard Model all these processes proceed through loop di- 
agrams and thus are relatively suppressed. In the exten- 
sions of the Standard Model the contributions stemming 
from the diagrams with " new" particles in the loops can 
be comparable or even larger than the contribution from 
the SM. Thus getting experimental information about 
rare decays puts strong constraints on the extensions of 
the SM or, if we are lucky, can lead to a strong disagree- 
ment with the SM predictions providing an evidence of 
some " new physics" . To make a rigorous comparison be- 
tween experiment and theory, one has to get refined the- 
oretical predictions for the rare decay at hand. In par- 
ticular, perturbative QCD corrections through next-to- 
next-to leading order (NNLO) in as are needed. When 
calculating higher order QCD corrections, along with the 
calculation of virtual corrections, it is necessary to take 
into account real emission of gluons and quark-antiquark 
pairs. Then the phase space integrals contain infrared 
and coUinear singularities, which can be regularized us- 
ing dimensional regularization scheme. 

Sometimes however it is more convenient to introduce 
a small mass as an infrared regulator which can regularize 
both the soft and coUinear divergences. This can be for 
instance massive s quark for 6 — > 57 [l| and h syj 

or massive light quarks for b — sqq-f {q = u,d,s) Q. 
Another example is the calculation of CP- asymmetry and 
width difference-Ar in B^^d — Bs,d mixing at NLO and 
NNLO in as , when the calculation gets simplified by the 
usage of a small gluon mass as an infrared regulator. 

In this article we derive formulae for the 1-^3 4] and 

1 — > 4 phase spaces, where the final state particles are 
massive. Following the approach of [5| for the 1 — )■ 3 
process we derive formulae for the phase space in the 
rest frame of two particles in the final state and obtain 
a factorized parameterization of the phase-space for the 
case when all three particles in final state are massive. 
In case of 1 — ;> 4 process we consider similar method, 
deriving the phase space formula in the rest frame of 
three particles in the final state. Here we come to the 
factorized parameterization of phase-space when only one 



particle in the final state is massive. 

The paper is organized as follows. In Section II we dis- 
cuss three-particle phase space. Section III is devoted to 
the investigation of four-particle phase space. In Section 
IV we apply our technique to the calculation of master 
integral connected with the B — B mixing at 0{a^). The 
Appendix contains formulae of 1 — > 3 and 1 — > 4 inte- 
grated phase spaces for some particular cases. 



II. THREE-PARTICLE PHASE SPACE 

Here we consider the three particle decay, when par- 
ticle with momentum p decays into three particles with 
momenta and masses m^, z = 1, 2, 3; p = P1+P2 +P3- 
We start from the well-known expression for the differ- 
ential decay width: 
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dT = IMI^D^. 
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where |Mp is the squared matrix element, summed and 
averaged over spins and colors of the particles in the fi- 
nal and initial states respectively and m is the mass of 
decaying particle. 

The phase space formula in the rest frame of particles 
1 and 3 with momenta Pi and p^ was derived in Q : 
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where d = 4— 2e and " dimensionless" momenta p = p /m, 
Pi = Pi/"^-, i=li2,3 are introduced. In (12.11) we introduced 
an additional integration over S13 = q^ , q — Pi + Ps 
(in the considered frame q = (y^iis, 0) ). Now we can 
integrate over d — 1 components of pi and p2 using the 
spatial parts of two d-dimensional S functions. To carry 
out the remaining integrations we choose the coordinate 
axes in a way that particle momenta have the following 
components: 



p = (£;,|;^,0,0,....), 

P3 = (£^3, |P3|cost9, |p3|sini?,0. 



(2.2) 
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where the dots correspond to the components of extra 
space dimensions, which are all zero. Making use of the 
remaining two one-dimensional i5-functions we express E 
and E3 in the following way: 



1 + 513 - X2 



X3 + 513 - Xi 



(2.3) 



where Xi 



II 2 



1,2,3 are dimensionless vari- 



ables. After integration over those angles on which |Mp 
does not depend, we obtain a factorized formula for the 
three massive particles phase space (unlike 4] , where the 
corresponding formula is not factorized): 



D$(l -> 3) 



^2(d-3) (d-4)/2 

24 (27r)2'^-3 "^13 

^2(d-3) (d-4)/2 



24 (27r)2d-3 13 
where fi^ is the solid angle in d dimensions. 



IpsI sini?)'' ^ d-d dsn 



^(l-cos2i?)"^dcosi?dsi3 



(2.4) 



V{d/2) ' 



(2.5) 



and 



b1 = 



1 



1 



iC2 - Si3j 



4X2 



2^/s 
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; v/(2;i + 2:3 - Sl3)^ - 4X1X3 (2.6) 



are 3-momenta expressed by Kdllen functions. 
For the integration limits we have 

(V^+ ^/^)^ < si3 < (1 - ^ f , 

-1 < C0S1? < 1. 



(2.7) 



We represent all scalar products via variables S13 and 
cosi^: 



P-P2 = 



1 + X2 - S 
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P-Pl = 



2(1 + Xi +X3)-X2-X 



P-P3 
P2-P1 = 
P2-P3 
P1-P3 



-2(xi-hx3)-X2 + 2si3 + X 



2 + 4x3 - X2 - 2si3 - X 



-X2 - 4x3 + X 



Sl3 - Xi - X3 



where 



X 



(2.8) 



-4|p1 |p3| cosz9 + 1 + 2(xi + X3) - si3 

(1 + Si3 - X2)(xi - X3) 



Sis 



(2.9) 



and IpI and \p3\ are defined in (2.6). 

From the formula (|2.4p for the three-particle phase 
space one can get a simpler formula assuming that the 
masses of two particles in the final state are equal to 
zero: 7712 — 7713 = 0. In that case we make the replace- 
ments: COST? = 2Ai — 1 and S13 = A2(l — xi) -f xi, where 
< Ai,2 < 1. Then we get: 



/^^(l ^ 3) 



2(d-3). 



f^d-2 
2d+l (27r)2'i-3 



dXi d\2 



(2.10) 



X (l-xi)^"-^ [A2(l-Xi)+Xi 



(2-d)/2 



[A2 (1 - A2)]'^-' [Ai(l-Ai)]^ 



If we put xi = (i.e. ttii = 0) in the expression (|2.10l) . 
we get the formula (2) of @. 



In the Appendix A we study the final state mass ef- 
fects in some interesting phenomenological cases, provid- 
ing fully integrated analytic formulae for the 1 — >■ 3 pro- 
cesses. 



III. FOUR-PARTICLE PHASE SPACE 



Now we proceed to the case of four particle decay. For 
this case one can derive the expression for the phase space 
in the similar way as in The only difference is that 
instead of two particle's rest frame introduced in here 
the rest frame of three particles is used (see Q). We 
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consider a particle with momentum p , which decays into 
four particles with momenta and masses 771^, i = 1, 4 
and p — Pi + P2 + P3 + Pa,- Here also we introduce 
" dimensionless" momentap = p /m, pi — Pi/m, i=l,...,4. 
In the rest frame of particles 2, 3 and 4 with momenta 
P2> ^3) Pi (where P2 + P3 + Pa = 0) four- particle phase 
space can be written in the following way 



D$2 = 



(S234) 



(d-2)/2 '± P_ 



Pi 



oS'ip- 



P2 



2p° {2nY-'^2p1 
d'^-^Ps, d'^-^pi 



(3.1) 

Pi - 9) , 



2'kY^^2pI {2TTf-^2pl (27r)'i-i2pO 

X5'^{q-P2-P3-PA), 

where S234=(?'2+P3+P4)'^ and g is a d-dimensional vector 
with components (-^5234, 0) in the rest frame of particles 
2, 3 and 4. We choose the coordinate axes in such a 
way, that the particle momenta have the following com- 

I 



ponents: 

P2 = (£;2,|P2|, 0,0, ....), 

Pi = (£;i,|pi|cosi?i,|pi|sin7?i,0, ....), (3.2) 

P3. = (£'3,|p3|cOSl?3,|p3|sini93 C0S(?!)3, |p3|sin7?3Sin03,....), 

where the dots, as before, correspond to the components 
of extra space dimensions, which are all zero. We can 
immediately integrate over those angles on which none 
of the scalar products depends and obtain 



d''-^P2^nd-l\p2\''-^ d\p2\, 

d'^-^pi = r!d_2 {&\ni)iY-^ dcos^i \pi\''-^ d\pi\ 
d''-^P3 = {smcl>3f-^ d cos 03 {sm^s)"-^ 

d-~2 



(3.3) 



xd cos 793 Wsr' d\p3\. 

Pa and p can be integrated out immediately using the 
spatial parts of the 6 functions in (3.1). Thus for the 
phase space we get 



where 



(s234)('^-')/2 nd-2{sm^iy-^ dc0S7?i jpil'^^^ 



5{E-El- , 



2E {2t:Y-'^2Ei 
^d-i\p2\''-^d\p2\ nd_3(sin03)''-'dcos(/)3(sin793)''^^dcos793|P3|''"'d|P3| 5 {V^-E2 -E-s- Ea] 



{2TTf-^2E2 



(27r)'^-i2£;3 



{2itY-^2Ea 



(3.4) 
-, (3.5) 



Ea = \/\p2^ + IpsP + 2|p2||p3| cos'(93 + xa , (3.6) function in d$i and (3.7) we get: 

(3.7) 



Xi = VTk^lvr? , i = 1, 4 



\-Xx- S234 p 1 - Xi + S234 
-^1 — , , — 



2VS234 



2^/5234 



(3.8) 



where = \J El — Xi , i = 1, 2, 3. Making use of the re- 
maining two one-dimensional 5-functions, we can express From the 8 function in d<J>2 together with (3.6) we find 
Ex and E-^ in terms of the other variables. Using the 5 the expression for En, 

I 



En 



1 



''2{yf^A~E2) +2|p-2PcOs2,?3 



{{E2 - VS234) (S234 - 2i;2Vs234 + X2 + X3 - Xa) 



+ IP2I COSl^sy (S234 - 2£;2\/s234 + 2:2 - 2^3 - 3^4) - 4x3 (x4 + |p2 P siu^ 193)) 



(3.9) 



Finally, introducing new variables: z\ = cosi?!, Z31 — cos7?3, Z32 = cos(/)3, we get for the four-particle phase space: 



£)$(1 4) 



32 



(i-3 



(2^) 



3d-4 



dz\ dzzi dz32 dE2 ds234 



X (1-Z?)-^(1-Z3\)^(l-Zf2)^ 



1-3 



r^A-E2 + i^bll ^31 



(3.10) 



where 



= \Pl\ \P2\ \P3\ \JS23A ■ 



The integration variables 3234,-^2,2^1,2^31,-232 in p.lOp 
(3 11) have the following limits 

(V^+ V^)' < S234 < (1- V^)', (3.12) 
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g234 + X2 - {^/xj + ^/xlY 

2^*234 
-1 < z, < 1 . 



(3.13) 



Different possible scalar products between 4-monienta 
of the particles can be obtained automatically from the 
expressions for their components (3.2) using the conser- 
vation of momentum. All the scalar products can be 
expressed through the 5 integration variables in (I3.10p : 



Pi.p2 ^ EiE2~ 
P2-P3 



|Pl||P2| C0S1?i , 
_ _ S2£^3-b^2||p3|cOS1?3, (3.14) 

Pi.p3 ^ EiEs- |pl||p5|(cosi?i C0SZ?3 
+ sin i?i sin cos 03 ) , 



In right-hand side of (3.15) the expressions given in (3.14) 
should be used. Note that in (3.14) and (3.15) Ei and 
i?3 are determined by formulae (3.8) and (3.9). 

We can also obtain a factorized formula for the 4 par- 
ticle phase space with 1 massive and 3 massless particles 
in the final state, taking X2 — — ~ 0. Making the 
following substitutions in p.lOp 



S234 = (1 - ^AT)'Al, E2 = I^/X'iil - A2)(l - V^) 



A2(l-A4)-A4 

^31 = -T-r-. , N , ^ , zi 



and 



2^32 



A2(l 
2A5- 



- A4) + A4 ' 
1, 



2A3 



(3.16) 



PI.P4 = Eiy/3234:- P1-P3 -~Pl-P2 
P2-Pi = E2^S23i - X2 - P2-P3 , 
PS.Pi = i?3%/S234 - X3 - P2.P3 , 



the differential phase space can be expressed in the fol- 
(3.15) lowing form: 



,3d-8 



D$(l 4) = 



dA 1 (iA2 d A3 (iA4 (iAs 



128 (27r)3^-4 
X (1 - V^)3'^-7 [(1 _ Ai)((l + V^)2 - Ai(l - V^)2)] ^^-^'^'^ 
X [Ai(l - A2)]'-' [A2(l ~ A3)A3(1 - A4)A4]''-')/' [A5(l - \,)f~'^'^ . 



(3.17) 



where A^ = [0,1], i = 1...5. The corresponding scalar S134 = (pi + P3 + ^4)^ = (sjj4 — Si34)A3 -I- S134, 

products in terms of the Ai will have the following form _ + Pzf = (s]*3 — •''j^3)A5 + s^g (3 18) 

S234 = {P2+P3+Pif ^ {'^- y/xifXi-, 

S34 = (p3 +^4)^ = (1 - \/2?r)^AiA2, 

S23 = (P2 +P3)^ = (1 - V^)^Ai(l - A2)A4, where 



'134 



1 + A2 - Ai(l - A2)(l - V^)2 + xi(l - A2) ± (1 - Vir)(l - \2)^{l - Ai)((l + ^/^)2 - Ai(l - V^)2) 



Si3 = a;i + - (1 - Vil) [(A2 (1 - A4) + A4) (1 + - Ai(l - ^)) 



+ (A2(l - A4) - A4) (1 - 2A3) ^(1 - Ai) ((1 + ViT)2 - Ai(l - V^)2) 
±2(1 - V^)y^A2(l - A3)A3(1 - A4)A4(1 - Ai) ((1 + V^)2 - Ai(l - V^)2) . 



(3.19) 



r 



As it is expected if we put xi = in p.l7p we get the In the Appendix B we consider the final state mass 
formula (21) of @ derived for massless particles in the effects in some interesting phenomenological cases, pro- 
final state. viding fully integrated analytic formulae for the 1^4 

processes. 
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FIG. 1: Left: an example of an infrared divergent diagram 
contributing to the B ~ B process at NNLO. The virtual glu- 
ons are represented by curly lines. Right: one of the master 
integrals appearing after reduction of the left diagram with 
two cuts A and B (c- is a charm, g- is a gluon and i- is either 
charm or beauty or a light quark). 



IV. CONTRIBUTION TO B 

0{al) 



B MIXING AT 



In this Section we illustrate how one can use formulae 
derived in previous sections for calculation of the absorp- 
tive part of a B — _B mixing at 0{al), which contributes 
to the width differences AF and the CP-asymmetry. 

An example of a diagram that contributes to the B — 
B mixing at ©(a^) is shown in Fig. 1. This diagram 
contains both 3-particle and 4-particle cuts. The light 
quark masses we put equal to zero. 

We introduce a small mass rUg for a gluon as an in- 
frared regulator. In this way, in the matching one doesn't 



need e and parts of NLO and LO Wilson coefficients 
correspondingly. Moreover, in the Effective Theory side 
the renormalization is much easier, because the (1/e") 
terms appear only as UV singularities and so their cance- 
lation can be tracked easier. Further in order to eliminate 
IR finite terms coming from the interference of IR (1/e) 
with UV structures of 0(e) a special class of operators 
(oc e) has to be introduced in case of dimensional regular- 
ization. Whereas with introduction of a gluon mass those 
terms vanish automatically. This is the major difference 
in using dimensional regularization or gluon mass for IR 
singularities. 

The reduction to the master integrals is done by means 
of a program FIRE 8] . 

Here we consider one of the most complicated master 
integrals (MI), shown on the right side of the Fig. 1. This 
MI has one three-particle cut and one four-particle cut. 
While in the case of 3-particle cut i on the right side of 
the Fig. 1 denotes either heavy rric, mf, or light quarks, in 
the case of a 4-particle cut i stands for only light quarks. 
The corresponding cuts can be calculated analytically as 
an expansion over a gluon mass. The coefficient of a 
corresponding MI contains inverse masses of nig up to 
the forth power. 

The 3-particle cut (cut A) can be factorized into a 
product of a qq fermion loop and a 3-particle phase space, 
which is presented in the Appendix, formula (jA.l[) . In the 
case of 4-particle cut we get: 



TB = ni [I67r2(2 - 3a;,)x, + \%2x\\Jx^(a) - 96(2 - x,)x,\Jy2 (a^) 

+ 12V1 - 4a;c(2a;c + 1)(8 log(l - <t) - 4 log(fT) - 5) 
-H48 log(a) {xl log(o-) - 4(2 - x^)xr_ log(l - cr) + (1 - x^f + 2(2 - Xc)x^ log(l - 4xe)) 
^Xg ■ (-87r2(5 - 6xc) + 96(1 - 2xc)U2{(7) + 48(3 - 2xc}U2 (cr^) 

-h241og(a)(4(3 - 2xc) log(l - a) + {1 - 2xc) log(cr) - 2(3 - 2xc) log(l - 4xc)) 
+8V1 - 4:Xc (77r2 - 6 - 48 log^(l - a) - 31og^(cr) + 481og(l - a) \og{a))) 
4 



(-28t:^xJ. 



3 6xt 



lOXr 



■log(a)) 



3 (l-4a;e)3/2 

-hl21og(f7) (1 - 2xc - 2x1 - VI - 4a;c(l - x^ 
-24 log(l - cr) (l - 2xc - 2x1 - log(l - (t) + 8a;^ log(cr))) 
- \og{xg) ■ (-24V1 - ^x^{2xc + 1) - 96a;c(l - Xc) log(a) 

+xg (96(1 - Xc) log(a) + 96V1 - 4ic (2 log(l - g) - log(a))) 



24 



3 (1 - 4Xc)3/2 



1^1 - 2x^ ~ 2x1 - 16a;2 log(l - cr) -f ^xl log(cr)) 



-I- log^(a;g) • -Xg24Vl - 4a;c + ; 



48x2 



(1 - 4xe)3/2 



+ 0{xl) 



(4.1) 



where n; is the number of light quarks (n; = 3), Xc = {rric/mi,)'^ , Xg ~ {rng/rri},)'^ and cr 
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In the expression (4.2) the infrared singularities ap- 
pear as log"(xg), n—1,2. The infrared singularities of 
the diagram shown on the left side of Fig. 1, together 
with the log"'(a;g)-s of the other diagrams get canceled 
in the matching with the correspondin g E ffective The- 
ory diagrams with Q, Qs insertions jol. Il0l|. It must be 
mentioned also, that the diagrams both in the full and 
Effective Theory have to be renormalized. In addition, in 
the Effective Theory corresponding Evanescent operators 
have to be taken into account as well [ill, [13 • 

V. CONCLUSIONS 

To conclude, we have presented new formulae for com- 
puting 1 — )■ 3 and 1 — >■ 4 processes, with massive particles 
in the final states. We have also obtained a factorized for- 
mula for the 4 particle phase space with 1 massive and 
3 massless particles in the final state. On the example 
oi B — B mixing we demonstrate the capability of our 
technique for NNLO calculations, considering one of the 
most complicated master integrals. Further we study the 
final state mass effects in some interesting phenomeno- 



logical cases providing fully integrated analytic formulae 
for 3 and 4 particle phase spaces. 
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APPENDIX A: FULLY INTEGRATED 
3-PARTICLE PHASE-SPACE 

Here, as an application of the formula \2A\ we present 
fully integrated 3-particle phase space. For the case of 
TOi = 7712, TOa = miight {mught < ^1,2), expanding ((2^ 
over miight I'm and e and performing integration, we get 
in the MS scheme 



PSsinii, nil, miight) 





-2£ 2 r 




771 


m 


1287r3 



-Vl - 4x1(1 + 2xi) + 2x1 (1 - xi) log(cr) 



+Xi (-2 (1 - xi) log(a) + VI - 4x1 (log(x£) - 4 log(l - a) + 2 log(a))) 

2{1 ~lxi)y^ (-1 + 2x1 + 2x1 - ^4 Mxi) - 41og(l - a) + 21og(a))) 



+ e ■ 



^(1 - xi)xi + 4(1 - xi)xi (Li2(CT) + Li2 (a^)) 

+(1 + 2xi)Vl-4xi ( ^ - ^ log(l - 4xi) - 2 log(l - a) 



log(CT; 



-(1 - xi) log(cr) (1 - 5xi -I- 6x1 log(l - 4xi) - 12xi log(l - cr) + Xi log(cr)) 
-Xi ( — (1 - xi) - 4(1 - xi) (Li2(a) + Lis {<7^)) 

+ (1 - xi) log(CT) (-6 + 6 log(l ~ 4xi) - 12 log(l -a) + \og{a)) 
+V1-4X1 (41og(l - c7)(21og(l - <j) - 21og(a) - 3) + (6 - log((7)) log(a) 



+xj 



+ log(l - 4xi)(41og(l - a) - 21og((T) - log{xe)) 

1 /. o o 1 



■log^(x^)-|-31og(x£)-27r2 + 2 



(1 -4xi)3/2 



47r2.x? 



(22x? - 14x1 + 7) 



1 



+ \og{xe) i^xl \og{xe) + - {-6x1 + ^xi - l) + 2x? log(l - 4xi) 
-2 log(cr) {-xj log(cr) + 2xj - 2x1 + 1 - Vl -4xi(l - xi)) 
+ i log(l - 4xi) {-16x1 log(l - cr) + 8x1 log(f^) - 2x? - 2xi + l) 
-Mog(l - a) (4 {2x1 - + 1) - 16x2log(l - a) + lex^ log((7)))]] +0{€^,xj) 



(A.l) 
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For the case mi > and m2 = we have 



PSsimi, 0,1713) 



1287r3 



^{1 + xi+ a;3)V'T^-^(2?7Ta?3yT|a?7^^X3p 



-2(xi + X3 — 2x1X3) coth 



^ 1- V^)' 



+ (•^1 -2:3) (^2 log 



ilog(16xiX3) ) 



For mi 3> m3 = mught, expanding over mught/rn, we get also the terms proportional to 



PS 3(1111, Q^mught) 



-2e 



m 



1287r3 



(1 -x?) +xilog(a;i) 



+xe. ((1 - xi) log(a:£) - 2(1 - xi) log(l - xi) ~ xi log(xi)) + x^\- 



1 — Xi 
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^1 - ^i) ^ + 2^1(6 -xi- log(a:i)) log(a:;i) 



-21og(l - xi) (1 -x^ + a;ilog(a;i)) - 2xi(Li2(l - .xi) - Li2(xi)) 
+Xi ((1 - xi) (61og'(l - xi) - 61og(l - xi) + 31og(a;£) + 2) 

-i(l-xi)(41og(l-2;i)log(x,)+log'(a;,)) -TT^ (l - ^) 

-xi f-21og(l - X,) - + 3) log(xi) + 2(Li2(l - X,) - x.Li^ix,)) 



+x; 



l+Xi 



31og(l - xi) 



1 — Xl 

where the following notations were introduced 



log(a:£) ^ 7^ Xl log(a;i) 
1 — Xl 



xe = {miight/m)'^ , 2:1,3 = (mi,3/m)^ 

I 



1 - ^/T^^I^ 



1 + VI - 4X1 



(A.2) 



(A.3) 



(A.4) 



APPENDIX B: FULLY INTEGRATED 
4-PARTICLE PHASE-SPACE 

As an application of the formula p.lOp we give fully 
integrated 4-particle phase space. For the case when 



mi = m 2, m3 = m;ig/it, m4 = {niught < mi, 2), ex- 
panding (|3.10p over mught/'m and e and performing in- 
tegration we get 



PS4{mi,Tni,miight,0) 



81927r5 



Att^xI 1 



. -VI - 4x1 (1 + 20x1 -I- 12x^) -I- 4x1 (1 + Xl - 2xf) log{a) 

-^4x2 log(cr)(21og(l - 4x1 ) - 41og(l - a) + log(cr)) - 16x?Li2(-(T) 
+xe ■ (-16x^Li2(cr) + 8x1(2 - xi)Li2 (cr^) + 4 (x? - 1 + VI -4xi(2xi + 1)) log(CT) 
-4xi log^(CT) - log(l - a) (8VI -4xi(2xi + 1) - 16(2 - xi)xi log(cr)) 

-8(2 - xi)xi log(CT) log(l - 4xi) + log(x£) (8(1 - xi)xi log(cr) + 2V1 -4xi(2xi + 1)) 

4 , 

--7r2(2 - 3xi)xi + 3Vl-4xi(2xi + 1) 

+xl ■ (-4(1 - Xl) log(a) - VI - 4x1 (3 + 8 log(l -a)- 41og(a) - 2 log(x,)))] + O (e, x^) 



(B.l) 
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For mi > TO3, 7712 = TO4 = we have 



P54(mi,0,m3,0) 



81927r5 



4 (xl - xl) tanh ^ 



+ i Vl - 2(a;i + 2:3) + {xi - 2:3)2 (1 + + a;^ + 10(a;i + X3 + xia;3)) 
+2{xi - X3) log(4a;i) + {xl{l - 2x3) + 2x3(1 - X1X3) + xj + Axix^ log(x3)) 
X {\og{AxiX3) - 2 log {1-X1-X3 + v/l - 2{xi + 2:3) + [xi - X3)2^ ^ 

-4(a;i - 3:3) log {l + xi~ xj, + \Jl- 2{xi + 2:3) + (xi - 2:3)2^ 

2xi \ 



— 8a:ia:3Li2 — 

^2:1+ 2:3 -1-^1- 2(a;i + 2:3) + [xi - X3) 

, „ -r • / a;i + 2:3 - 1 - ^1 - 2(2:1 + X3) + {xi - xz) 

+82;iX3Li2 

\ 22:3 



For mi ^ TO3 = niiigfit, expanding the expression (|B.2p over m3, we get 



PS4{mi,0,miight,0) = 



81927r5 



- a;i) {xl + lOxi + 1) + 22;i(.Ti + 1) log(2;i) 

+Xi ((1 - xl) (3 - 41og(l - 2:1) + 21og(2;f)) + 22;i log(2;i)(-2;i - 41og(l - 2:1) + 21o^ 
-8a:iLi2(l-2;i)) 

+2;2(-(l - 2;i)(3 + 41og(l - 2:1) - 2 log(a:,)) - 2x1 log(2;i))] + O (e, x^) . 

I 



[xi)) 



(B.2) 



(B.3) 
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